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Abstract: The packet is the fundamental unit of transportation in modern communication 
networks such as the Internet. Physical layer scheduling decisions are made at the level of 
packets, and packet-level models with exogenous arrival processes have long been employed to 
study network performance, as well as design scheduling policies that more efficiently utilize 
network resources. On the other hand, a user of the network is more concerned with end-to-end 
bandwidth, which is allocated through congestion control policies such as TCP. Utility-based 
' flow-level models have played an important role in understanding congestion control protocols. 

> ^ , In summary, these two classes of models have provided separate insights for flow-level and 

' packet-level dynamics of a network. 

In this paper, we wish to study these two dynamics together. We propose a joint flow-level 
and packet-level stochastic model for the dynamics of a network, and an associated policy 
for congestion control and packet scheduling that is based on a-weighted policies from the 
literature. We provide a fluid analysis for the model that establishes the throughput optimality 
of the proposed policy, thus validating prior insights based on separate packet-level and flow- 
level models. By analyzing a critically scaled fluid model under the proposed policy, we provide 
CZ2 , constant factor performance bounds on the delay performance and characterize the invariant 

^ ' states of the system. 

Keywords and phrases: Flow-level model, Packet-level model. Congestion control, Schedul- 
^ ' ing, Utility maximization. Back-pressure maximum weight. 
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1. Introduction 



' The optimal control of a modern, packet-switched data network can be considered from two distinct 

' vantage points. From the first point of view, the atomic unit of the network is the packet. In 

a packet-level model, the limited resources of a network are allocated via the decisions on the 
scheduling of packets. Scheduling policies for packet-based networks have been studied across a 
long line of literature (e.g., [30, 27, 24]). The insights from this literature have enabled the design of 
5^ I scheduling policies that allow for the efficient utilization of the resources of a network, in the sense 

of maximizing the throughput of packets across the network, while minimizing the delay incurred 
by packets, or, equivalently, the size of the buffers needed to queue packets in the network. 

Packet-level models accurately describe the mechanics of a network at a low level. However, 
they model the arrival of new packets to the network exogenously. In reality, the arrival of new 
packets is also under the control of the network designer, via rate allocation or congestion control 
decisions. Moreover, while efficient utilization of network resources is a reasonable objective, a 
network designer may also be concerned with the satisfaction of end users of the network. Such 
objectives cannot directly be addressed in a packet-level model. 
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Flow- level models (cf. [12, 3]) provide a different point of view by considering the network at 
a higher level of abstraction or, alternatively, over a longer time horizon. In a flow-level model, 
the atomic unit of the network is a flow, or user, who wishes to transmit data from a source to a 
destination. Resource allocation decisions are made via the allocation of a transmission rate to each 
flow. Each flow generates utility as a function of its rate allocation, and rate allocation decisions 
may be made so as to maximize a global utility function. In this way, a network designer can address 
end user concerns such as fairness. 

Flow-level models typically make two simplifying assumptions. The first assumption is that, 
as the number of flows evolves stochastically over time, the rates allocated to flows are updated 
instantaneously. The rate allocation decision for a particular flow is made in a manner that requires 
immediate knowledge of the demands of other flows for the limited transmission resources along the 
flow's entire path. This assumption, referred to as time-scale separation, is based on the idea that 
flows arrive and depart according to much slower processes than the mechanisms of the rate control 
algorithm. The second assumption is that, once the rate allocation decision is made, each flow can 
transmit data instantaneously across the network at its given rate. In reality, each flow generates 
discrete packets, and these packets must travel through queues to traverse the network. Moreover, 
the packet scheduling decisions within the network must be made in a manner that is consistent 
with and can sustain the transmission rates allocated to each flow, and the induced packet arrival 
process must not result in the inefficient allocation of low level network resources. 

In this paper, our goal is to develop a stochastic model that jointly captures the packet-level and 
flow-level dynamics of a network, without any assumption of time-scale separation. The contribu- 
tions of this paper are as follows: 

1. We present a joint model where the dynamic evolution of flows and packets is simultaneous. 
In our model, it is possible to simultaneously seek efficient allocation of low level network 
resources (buffers) while maximizing the high-level metric of end-user utility. 

2. For our network model, we propose packet scheduling and rate allocation policies where 
decisions are made via myopic algorithms that combine the distinct insights of prior packet- 
and flow-level models. Packets are scheduled according to a maximum weight policy. The 
rate allocation decisions are completely local and distributed. Further, in long term (i.e., 
under fluid scaling), the rate control policy exhibits the behavior of a primal algorithm for 
an appropriate utility maximization problem. 

3. We provide a fluid analysis of the joint packet- and flow-level model. This analysis allows 
us to establish stability of the joint model and the throughput optimality of our proposed 
control policy. 

4. We establish, using a fluid model under critical loading, a performance bound on our control 
policy under the metric of minimizing the outstanding number of packets and flows in the 
network (or, in other words, minimizing delay). We demonstrate that, for a class of balanced 
networks, our control policy performs to within a constant factor of any other control policy. 

5. Under critical loading, we characterize the invariant manifold of the fluid model of our con- 
trol policy, as well as establishing convergence to this manifold starting from any initial 
state. These results, along with the method of Bramson [4], lead to the characterization of 
multiplicative state space collapse under heavy traffic scaling. Further, we establish that the 
invariant states of the fluid model are asymptotically optimal under a limiting control policy. 
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In summary, our work provides a joint dynamic flow- and packet-level model that captures the 
microscopic (packet) and macroscopic (fluid, flow) behavior of large packet-based communications 
network faithfully. The performance analysis of our rate control and scheduling algorithm suggests 
that the separate insights obtained for dynamic flow-level models [12, 3] and for packet-level models 
[30, 27, 24] indeed continue to hold in the combined model. 

The balance of the paper is organized as follows. In Section 1.1, we survey the related literature on 
flow- and packet-level models. In Section 2, we introduce our network model. Our network control 
policy, which combines features of maximum weight scheduling and utility-based rate allocation is 
described in Section 3. A fluid model is derived in Section 4. Stability (or, throughput optimality) of 
the network control policy is established in Section 5. The critically scaled fluid model is described 
in Section 6. In Section 7, we provide performance guarantees for balanced networks. The invariant 
states of the critically scaled fluid model are described in Section 8. Finally, in Section 9, we 
conclude. 

1.1. Literature Review 

The literature on scheduling in packet-level networks begins with Tassiulas and Ephremides [30], 
who proposed a class of 'maximum weight' (MW) or 'back-pressure' policies. Such policies assign 
a weight to every schedule, which is computed by summing the number of packets queued at links 
that the schedule will serve. At each instant of time, the schedule with the maximum weight will be 
selected. Tassiulas and Ephremides [30] establish that, in the context of multi-hop wireless networks, 
MW is throughput optimal. That is, the stability region of MW contains the stability region of 
any other scheduling algorithm. This work was subsequently extended to a much broader class of 
queueing networks by others (e.g., [18, 7, 29, 6]). 

By allowing for a broader class of weight functions, the MW algorithm can be generalized to the 
family of so-called MW-a scheduling algorithms. These algorithms are parameterized by a scalar 
a G (0, oo). MW-a can be shown to inherit the throughput optimality of MW [14, 22] for all values 
of a G (0,oo). However, it has been observed experimentally that the average queue length (or, 
'delay') under MW-a decreases as a — )■ O"*" [14]. Certain delay properties of this class of algorithms 
have been subsequently established under a heavy traffic scaling [27, 6, 24]. 

Flow-level models have received significant recent attention in the literature, beginning with the 
work of Kelly, Maulloo, and Tan [12]. This work developed rate-control algorithms as decentralized 
solutions to a deterministic utility maximization problem. This optimization problem seeks to 
maximize the utility generated by a rate allocation, subject to capacity constraints that define a 
set of feasible rates. This work was subsequently generalized to settings where flows stochastically 
depart and arrive [17, 8, 3], addressing the question of the stability of the resulting control policies. 
Fluid and diffusion approximations of the resulting systems have been subsequently developed 
[13, 11, 32]. Under these stochastic models, flows are assumed to be allocated rate as per the 
optimal solution of the utility maximization problem instantaneously. Essentially, this time-scale 
separation assumption captures the intuition that the dynamics of the arrivals and departures of 
flows happens on a much slower time-scale than the dynamics of rate control algorithm. 

In reality, flow arrivals/departures and rate control happen on the same time-scale. Various 
authors have considered this issue, in the context of understanding the stability of the stochastic 
flow level models without the time-scale separation assumption [15, 10, 28, 20, 26]. Lin, Schroff, 
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and Srikant [15] assume a stochastic model of flow arrivals and departures as well as the operation 
of a primal-dual algorithm for rate allocation. However, there are no packet dynamics present. 
Other work [10, 28, 20] has assumed that rate control for each type of flow is a function of a local 
Lagrange multiplier; and a separate Lagrange multiplier is associated with each link in the network. 
These multipliers are updated using a maximum weight-type policy. In this line of work, Lagrange 
multipliers are interpreted as queue lengths, but there are no actual packet-level dynamics present. 
Further, these models lack flow-level dynamics as well. Thus, while overall this collection of work is 
closest to the results of this paper, it stops short of offering a complete characterization of a joint 
flow- and packet-level dynamic model. 

Finally, we take note of recent work by Walton [31], which presents a simple but insightful model 
for joint flow- and packet-level dynamics. In this model, each source generates packets by reacting 
to the acknowledgements from its destination, and at each time instant, each source has at most 
one packet in flight. Under a many-source scaling for a specific network topology, it is shown that 
the network operates with rate allocation as per the proportional fair criteria. This work provides 
important intuition about the relationship between utility maximization and the rate allocation 
resulting from the packet-level dynamics in a large network. However, it is far from providing a 
comprehensive joint flow- and packet-level dynamic model as well as efficient control policy. 

2. Network Model 

In this section, we introduce our network model. This model captures both the flow-level and the 
packet-level aspects of a network, and will allow us to study the interplay between the dynamics 
at these two levels. In a nutshell, flows of various types arrive according to an exogenous process 
and seek to transmit some amount of data through the network. As in the standard congestion 
control algorithm, TCP, the flows generate packets at their ingress to the network. The packets 
travel to their respective destinations along links in the network, queueing in buffers at intermediate 
locations. As each packet travels along its route, it is subject to physical layer constraints, such 
as medium access constraints, switching constraints, or constraints due to limited link capacity. A 
flow departs once all of its packets have been sent. 

In this section, we focus on the mechanics of the network that are independent of the network 
control policy. In Section 3, we will propose a specific network control policy to be applied in the 
context of this model. 

2.1. Network Structure 

Consider a network consisting of a set V of destination nodes, a set C of links, and a set of flow 
types. Each flow type is identifled by a fixed given route starting at the source link s{f) € C and 
ending at the destination node d{f) € V. At a given time, multiple fiows of a given type exist in 
the network, each flow injects packets into the network. 

The network maintains buffers for packets that are in transit across the network. At each link, 
there is a separate queue for the packets corresponding to each possible destination. Let £ = C xV 
denote the set of all such queues, with each e = {i, v) being the queue at link i for flnal destination 
V. Traffic in each queue is transmitted to the next hop along the route to the destination, and 
leaves the network when it reaches the destination. We define the routing matrix R € {0, l}^'^^ 
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by setting iige' — 1 if the next hop for queue e is queue e', and R^e' — otherwise. Traffic for a 
flow of type / enters the network in the ingress queue i{f) = (^s{f),d{f)^ € £. Define the ingress 
matrix T € {0, 1}^^-^ by setting T^f = 1 if L{f) = e, and Tg/ = otherwise. We will assume that 
the routes are acyclic. In this case, we can define the matrix 

E^(^I-R^y^ =I + R^ + [R^f + ... , (1) 

Under the acyclic routing assumption, Hg/g = 1 if and only if a packet arriving at queue e will 
subsequently eventually pass through queue e'. 



2.2. Dynamics: Flow-Level 



In this section, we will describe in detail the stochastic model for dynamics of flows in the network. 
The system evolves in continuous time, with t € [0, oo) denoting time, starting at t = 0. For each 
flow type / G J-", let Nf{t) denote the number of flows of type / active at time t. Flows of type / 
arrive according to an independent Poisson process of rate i^f. Flows of type / receive an aggregate 
rate of service Xf{t) G [0, C] at time t. Here, C > is the maximal rate of service that can be 
provided to any flow type. The total rate of service Xf[t) is divided equally amongst the Nf{t) 
flows. As flows are serviced, packets are generated. The evolution of packets and flows proceeds 
according to the following: 

• Packets are generated by all the flows of type /, in aggregate, as a time varying Poisson 
process of rate Xf{t) at time t. If Nf{t) = 0, then we require that Xf{t) = 0. 

• When a packet is generated by a flow of type /, it joins the ingress queue i{f) € 8. 

• When a packet is generated by a flow of type /, the flow departs from the network with a 
probability of < /ij < 1, independent of everything else. 

Thus, each flow of type / generates a number of packets that is distributed according to an inde- 
pendent geometric random variable with mean^ f/^/i ^^"^ ff^^ f^o^ departure process for flows of 
type / is a Poisson process of rate HfXf{t) at time t. We can summarize the flow count process 
Nf{-) by the transitions 

AT (t) ^ /^^^*) + ^ at rate z^/, 
■'^ \iV/(t) - 1 at rate fifXf{t). 

Define the offered load vector p € M;^ by pf = Uf/pf, for each flow type /. Without loss of 



generality, we will make the following assumptions:^ 

• p > 0, i.e., we restrict attention to flows with a non-trivial loading. 

• p < CI, i.e., we assume that the maximal service rate C is sufficient for the load generated 
by any single flow type. 

• HFp > 0, i.e., we restrict attention to queues with a non-trivial loadings. 



^The assumption that /i/ < 1 is equivalent to requiring that > 1, i.e., each flow is expected to generate more 

than one packet. This is reasonable since we require flows to arrive with at least one packet and for there to be some 
variability in the number of packets associated with a flow. 

^In what follows, inequalities between vectors are to be interpreted component- wise. The vector (resp., 1) is the 
vector where every component is (resp., 1), and whose dimension should be inferred from the context. 
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Denote by Aj(t) the cumulative number of flows of type / that have arrived in the time interval 
[0,t]. Denote the cumulative number of packets generated by flows of type / in the time interval 
[0,t] by Af{t). We suggest that the reader take note of difference between Aj(-) and Af{-). Let 
Df[t) denote the cumulative number of flows of type / that have departed in the time interval 
[0,t]. The evolution of the flow count for flow type / over time can be written as 

Nf{t) = NJ{S))+^f{t)-DJ{t). (2) 
2.3. Dynamics: Packet-Level 

As we have just described, flows generate packets which are injected into the network. These packets 
must traverse the links of the network from source to destination. In this section, we describe the 
dynamics of packets in the network. 

We assume that each queue in the network is capable of transmitting at most 1 data packet per 
unit time. However, the collection of queues that can simultaneously transmit is restricted by a set 
of scheduling constraints. These scheduling constraints are meant to capture any limitations of the 
network due to scarce resources (e.g., limited wireless bandwidth, limited link capacity, switching 
constraints, etc.). 

Formally, the scheduling constraints are described by the set S C {0, 1}^. Under a permissible 
schedule vr € 5, a packet will be transmitted from a queue e ^ 8 \i and only if vTe = 1. We assume 
that € 5. We require that each queue e be served by some schedule, i.e., there exists a vr G 5 with 
TTe = 1. Further, we assume that S is monotone: if a G 5 and a' € {0, 1}^ is such that a'^, < Ce for 
every queue e, then a' € S. Finally, denote by IT € {0, ij^^-^ the matrix with columns consisting 
of the elements of S. 

We assume that the scheduling of packets happens at every integer time. At a time r G Z-|_, let 
7r(r) G S denote the scheduled queues for the time interval [r, r + 1). For each queue e, denote by 
Qe(r~) the length of the queue e immediately prior to the time r (i.e., before scheduling happens). 
The queue length evolves, for times t £ [r, r + 1) according to^ 

Qe{t) = Qe{T~) - Vre(r)I{Q^(,-)>o} + ^^fi^f^^) " ^/(^")) + ^e'eVTe' (r)I{Q^, )>o} • 

Here, for each flow type /, Af{T^) is the cumulative number of packets generated by flows of type 
/ in the time interval [0,r). The term 7re(r)I|Q^(T-)>o} enforces an idling constraint, i.e., if queue e 
is scheduled but empty, no packet departs. Note that, over a time interval [r, r + 1), we assume the 
transmission of packets already present in the network occurs instantly at time r, while the arrival 
of new packets to the network occurs continuously throughout the entire time interval. 

Finally, let S-,^{t) denote the cumulative number of time slots during which the schedule tt was 
employed up to and including time r. Let Ze{T) denote the cumulative idling time for queue e up 
to and including time r. That is, 

T 

Ze{T) = ^e('S'7r(s) " S.^{s - 1)) I|q^(,)=o} • 

denotes the indicator function. 
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Then, the overall queue length evolution can be written in vector form as 



Q(t + 1) = Q(0) - (/ - i?^)n5(r) + (I - R^)Z{t) + VA{t + 1) 



(3) 



where we define the vectors 



Q{t) 4 

S{r) 4 [S.(t)] 



A{t) 4 [Af{t)\ 



ee£" 



3. MWUM Control Policy 

A network control policy is a rule that, at each point in time, provides two types of decisions: (a) 
the rate of service provided to each flow, and (b) the scheduling of packets subject to the physical 
constraints in the network. In Section 2, we described the stochastic evolution of flows and packets 
in the network, taking as given the network control policy. In this section, we describe a control 
policy called the maximum weight utility maximization-a (MWUM-a) policy. MWUM-a takes as a 
parameter a scalar a G (0, oo) \ {1}. 

The MWUM-a policy is myopic and based only on local information. Specifically, a flow generates 
packets at rate that is based on the queue length at its ingress, and the scheduling of packets is 
decided as a function of the effected queue lengths. 

At the flow-level, rate allocation decisions are made according to a per flow utility maximization 
problem. Each flow chooses a rate so as to myopically maximize its utility as a function of rate 
consumption, subject to a linear penalty (or, 'price') for consuming limited network resources. As 
in the case of a-fair rate allocation, the utility function is assumed to have a constant relative risk 
aversion of a. The price charged is a function of the number of packets queued at the ingress queue 
associated with the flow, raised to the a power. 

At the packet-level, packets are scheduled according to a maximum weight-a scheduling algo- 
rithm. In particular, each queue is assigned a weight equal to the number of queued packets to the 
a power, and a schedule is picked which maximizes the total weight of all scheduled queues. 

3.1. Control: Rate Allocation 

The first control decision we shall consider is that of rate allocation, or, the determination of the 
aggregate rate of service Xf{t), at time t, for flows of type /. We will assume our network is governed 
by a variant of an a-fair rate allocation policy. This is as follows: 

Assume that each flow of type / is allocated a rate Yf{t) > at time t by maximizing a (per 
flow) utility function that depends on the allocated rate, subject to a linear penalty (or, cost) for 
consuming resources from the limited capacity of the network. In particular, we will assume a utility 
function given a rate allocation of y > to an individual flow of type / of the form 



for some a E (0, oo) \ {!}. This utility function is popularly known as a-fair in the congestion 
control literature [19], and has a constant relative risk aversion of a. The individual flow will be 



A y 



l-a 



Vf{y) 



l-a 
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assigned capacity according to 



Yf{t) G argmax Vf{y) - Q'^(^f^{t)y. 



Here, represents a 'price' or 'congestion signal'. Intuitively, a flow reacts to the congestion 

in the network (or lack thereof) through the length of the ingress or 'first-hop' queue L{f). Then, 
if Nf{t) > 0, the aggregate rate Xf(t) allocated to all flows of type / at time t is determined 
according to 

Xj{t) = Njit)Yf{t) = argmax ——I Q'^^j^{t) x. 



x>0 



If Nf{t) = 0, we require that Xf{t) = 0. Further, we will constrain the overall rate allocated to 
flows of type / by the constant C . Thus, rate allocation is determined by the equation 



otherwise 



argmax — ^ Q^, (t) x if iV/ (t) > 0, 

x6[0,C] 



Given the strictly concave nature of the objective in this optimization program, it is clear that the 
maximizer is unique and Xf{t) is well-defined. 

Denote by Xf{t) the cumulative rate allocation to flows of type / in the time interval [0,t], i.e., 

Xf{t)^ fxf{s)ds. 
Jo 

Xf{-) is Lipschitz continuous and differentiable, since Xf{-) is always bounded by C. 
3.2. Control: Scheduling 

The second control decision that must be specified is that of scheduling. We will assume the fol- 
lowing variation of the 'maximum weight' or 'back-pressure' policies introduced by Tassiulas and 
Ephremides [30]. 

At the beginning of each discrete-time slot r G Z_|., a schedule tt{t) G 5 is chosen according to 
the optimization problem 



7r(r) G argmax VTe 



^XT-)-yRee'Q:^{T-) 



e'&S 



(4) 



argmax tt ' {I - R)Q''{t' 

7rG5 



where Q°{t~) = [Qe (''"")] gg^ ^ vector of component-wise powers of queue lengths, immediately 
prior to time r. In other words, the schedule '7r(r) is chosen so as to maximize the summation of 
weights of queues served, where weight of a queue e G £" is given by [{I — R)Q'^ {t~)] ^. Given that S 
is monotone, there exists a vr G 5 that maximizes this weight and is such that vrg = if Qe(7" ) = 0. 
We will restrict our attention to such schedules only. From this, it follows that the objective value 
of the optimization program in (4) is always non-negative. 

By the discussion above, it is clear that the following invariants are satisfied: 
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1. For any schedule vr and time r, St^{t) = St^{t — 1) if 

TT^il - R)Q"{t-) < a^{I - R)Q''{t'), 

for some cr € 5. 

2. For any queue e and time r, Z^ir) = 0. In other words, there is no idhng. 
4. Fluid Model 

In this section, we introduce the fluid model of the our system. As we shall see, the allocation of rate 
to flows in the fluid model resembles rate allocation of 'flow-level' models that has been popular 
in the literature [17, 3]. In that sense, our model on original time-scale operates at a packet-level 
granularity and on the fluid-scale operates at a flow-level granularity. 

4.1. Fluid Scaling and Fluid Model Equations 

In order to introduce the fluid model of our network, we will consider the scaled version of the 
system. To this end, denote the overall system state at a time t > by 

Z{t) ^ (^Q{t),Z{[t\),N{t),S{[t\),X{t),A{t),D{t),A{t)). 

Here, the components of the state 2{t) are the primitives introduced in Sections 2.2 and 2.3. That 
is, at times t € M+, we have. 



Q{t) 


_A 




where Qe{t) 


is 


the length of queue e. 


Nit) 


A_ 




where Nf[t) 


is 


the number of flows of type /, 


X{t) 


_A 




where Xf{t) 


is 


the cumulative rate allocated to flow type /, 


A(t) 


_A 




where Aj(t) 


is 


the cumulative arrival count of flow type /, 


m 


_A 




where Df[t) 


is 


the cumulative departure count of flow type /, 


A{t) 


_A 




where Af{t) 


is 


the cumulative packet arrival count of flow type / 



and, at times r € Z+, we have 



Z{t) = [Ze(T)]gg^, where Ze{T) is the cumulative idleness for queue e, 
S{t) = ['S'ttI''")] ^g^, where St^{t) is the cumulative time schedule tt is employed. 
Given a scaling parameter r € M, r > 1, define the scaled system state as 

^ (QW(t),zW(t),A^W(t),5(")(t),xW(t),A('^)(t),DW(t),AW(t)). (5) 

Here, the scaled components are defined as 

X(''){t)^r~'^X{rt), AW(t) ^r-^A{rt), 
A r-^D{rt), = r-^A{rt), 

Zir) (i) A ^-1 [(^i _ LrtJ )Z{ \rt] ) + ( \rt-] - rt)Z{ [rt\ )] , 
sir) A ^-1 [(^^ _ L^r^J )5( p^^^ ) ^ ( p^i^ _ L^^j )] _ 
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In the above, the components Z(-) and S{-) are hnearly mterpolated for technical convenience only. 

Our interest is in understanding the behavior of Z^^'\-) as r — t- oo. Roughly speaking, in this 
limiting system the trajectories will satisfy certain deterministic equations, called fluid model equa- 
tions. Solutions to these equations, which are defined below, will be denoted as fluid model solutions. 
The formal result is stated in Theorem 1. 

Definition 1 (Fluid Model Solution). Given fixed initial conditions q{0) G and n(0) G M;^, 
for every time horizon T > 0, let FMS(T) denote the set of all trajectories 

^ {q{t),z{t),n{t),s{t),x{t),ait),d{t),a{t)) 

G 3 = X X M;^ X M:^ X M;^ X M;^ X M;^ X M;^, (6) 

over the time interval [0, T] such that: 

(Fl) All components of i{t) are uniformly Lipschitz continuous and thus differentiahle for almost 

every t G (0, T) . Such values of t are known as regular points. 
(F2) For all t G [0,T], n{t) = n(0) + a(t) - d{t). 
(F3) For all t G [0,r], a(t) = vt. 
(F4) For all t G [0,r], = diag(/x)x(t). 
(F5) For all t G [0,T], 

q{t) = q{0) -{I- R^)Us{t) R^)z{t) + Ta{t). 

(F6) For all t G [0,r], a{t) =x{t). 
(F7) For all t G [0,r], 1^ s{t) = t. 

(F8) Each component of z{-), s{-), and x{-) is non- decreasing. 

In addition, define the set FMS°(T) to be the subset of trajectories in FMS(T) that also satisfy: 
(F9) If t ^ (0, is a regular point, then for all f , 

{x^-°'n°i(t) 
argmax — ^ gl^^^^ (t) x if 7i/ (t)>{), 
Vfl\Xj {— pf) otherwise, 

where^ Xf{t) = Xf{t). 
(FIO) If t ^ (0, r) is a regular point, then for all vr G 5, ^^^(t) = 0, if 

TT^{I - R)q''{t) < max {I - R)q''{t). 

(Tg<S 

(Fll) If t G (0, r) is a regular point and nf{t) = for some f G T, then q^(^f^{t) = 0. 
(F12) For all t G [0, T], z{t) = 0. 

Note that (F1)-(F8) correspond to fluid model equations that must be satisfied under any 
scheduling policy, and, hence, are algorithm independent fluid model equations. On the other hand, 
(F9)-(F12) are particular to networks controlled under the MWUM-a policy. (F9) captures the 
long-term effect of the rate allocation mechanism through the a-fair utility maximization based 

"We use the notation e{t) to denote f^e{t) for 6: [0, T] R. 
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policy. Indeed, in a static resource allocation model, (F9) can be thought of as the primal update 
in an algorithm that seeks to allocate rates to maximize the net a-fair utility of the flows subject 
to capacity constraints. (FIO) captures the effect of short-term packet-level behavior induced by 
the scheduling algorithm. Specifically, the characteristics of the MW-a packet scheduling algorithm 
are captured by this equation. 

4.2. Formal Statement 

We wish to establish fluid model solutions as limit of the scaled system state process Z^^\-) as 
r — > oo. To this end, fix a time horizon T > 0. Let D[0,T] denote the space of all functions from 
[0,r] to 3j as defined in (6), that are right continuous with left limits (RCLL). We will denote 
the Skorohod metric on this space by ci(-,-) (see Appendix B for details). Given a fixed scaling 
parameter r, consider the scaled system dynamics over interval [0, T]. Each sample path Z^'^\-) of 
the system state is RCLL, and hence is contained in the space D[0, T]. 

The following theorem formally establishes the convergence of the scaled system process to a 
fluid model solution of the form specified in Definition 1. 

Theorem 1. Given a fixed time horizon T > 0, consider a collection of scaled system state processes 
{Z^'''\-) : r > 1} C D[0, T] under an arbitrary control policy. Suppose the initial conditions 

lim Q('^)(0) = q{0), lim iVW(O) = n(0), (7) 

r— >oo r—>oo 

are satisfied with probability 1. Then, for any e > 0, 

liminf P(z^''\-) e FMSe(r)') = L 

r->oo \ / 

Here, FMS£{T) is an e-flattening of the set FMS(T) of fluid model solutions, i.e., 

FMS,(T) ^ {x e D[0,r] : d(x,y) < y G FMS(r)}. 
Additionally, under the MWUM-a control policy, we have that 

liminf pfzW(-) e FMS"(r)') = 1. 

Here, FMS"(r) is an e-flattening of the set FMS"(r) of MWUM-a fluid model solutions, i.e., 

FMS-(r)^{xeD[0,r] : d(x,y)<e, y G FMS°(r)}. 

Theorem 1 can be established by following a somewhat standard sequence of arguments (cf. 
[4, 25, 13]). First, the collection of measures corresponding to the collection of random processes 
{Z^'^\-) : r > 1} is shown to be tight. This establishes that limit points must exist. Next, 
it is established that each limit point must satisfy the conditions of a fluid model solution with 
probability 1. The tightness argument uses concentration properties of Poisson process along with 
the Lipschitz property of queue length process. A detailed argument is required to establish that, 
with probability 1, the conditions of fluid model solution are satisfied. The complete proof of 
Theorem 1 is presented in Appendix B. 
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5. System Stability 

In this section, we characterize the stabihty of a network under the MWUM-a pohcy. In particular, 
we shall see that the network evolves according to a Markov process is positive recurrent as long 
as the system is underloaded. In other words, the system is maximally stable. In order to construct 
the stability region under the MWUM-a policy, first define the set of packet arrival rates by 

A = |a € Mf. : 3 s € Rf with A < Us, l^s < l} . (8) 

Imagine that the network has no packet arrivals from flows, but instead has packets arriving ac- 
cording to exogenous processes. Suppose that A G is the vector of exogenous arrival rates, so 
that packets arrived to each queue e at rate Ag. Then, it is not difficult to see that the network 
would not be stable under any scheduling policy if A ^ A. This is because there is at least one 
queue in the network that is loaded beyond its service capacity. Hence, the set A represents the 
raw scheduling capacity of the network. 

The set A can alternatively be described as follows: Given a vector A € M^, consider the linear 
program 

PRIMAL(A) = minimize l^s 

s 

subject to A < lis, 
s e Mf . 

Clearly A G A if and only if PRIMAL(A) < 1. The quantity PRIMAL(A) is called the effective load 
of a system with exogenous packet arrivals at rate A. 

Now, in our model, packets arrive to the network not through an exogenous process, but rather, 
they are generated by flows. As discussed in Section 2.2, each flow type f € T generates packets at 
according to an offered load of The generated packets are injected into the network according 
to the ingress matrix F, and subsequently travel through the network along pre-determined paths 
specified by the routing matrix R. Let A S be the vector of implied loads on the scheduling 
network due to the packets generated by flows. It seems reasonable to relate A and the vector 
p € M;^ of offered loads according to A = Fp + X. Equivalently, we define A = HFp, where H is 
from (1). We define the effective load L(p) of our network by L{p) = PRIMAL(Srp). 

Given the above discussion, it seems natural to suspect that the network's scheduling capacity 
allows it to operate effectively as long as L(p) < 1. This motivates the following definition: 

Definition 2 (Admissibility). A vector p G M;^ of offered loads admissible if\-{p) < 1. Similarly, 
p is strictly admissible if \-{p) < 1. Finally, p is critically admissible if \-{p) = 1. 

We will establish system stability, or, more formally, positive recurrence, when offered load is 
strictly admissible. To this end, recall that the system is completely described by the Z{-) process. 
Under the MWUM-a policy, the evolution of all the components of Z{-) is entirely determined by 
(A^(-), Q{-)) . Further, the changes in (A^(-), Q{-)) occur at times specified by the arrivals of a (time- 
varying) Poisson process. Therefore, tuple (A^(-), (5(')) forms the state space of a continuous-time 
Markov chain. The following is the main result of this section: 

Theorem 2. Consider a network system with strictly admissible p operating under the MWUM-a 
policy. Then, the Markov chain {^N{-),Q{-)^ is positive recurrent. 
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It is worth noting that if \-{p) > 1, then at least one of the queues in the network must be, 
on average, loaded beyond its capacity. Hence the network Markov process can not be positive 
recurrent or stable. 

The proof of Theorem 2, provided next in Section 5.1, uses a fluid model approach. Dai [5] 
pioneered such an approach for a class of queueing networks. However, this result does not apply 
to the present setting, and a specialized analysis is needed. Conceptually, the fluid model approach 
involves two steps: (1) derive strong stability of fluid model, and (2) use strong stability to establish 
the positive recurrence of the original Markov chain. To this end, define Lyapunov function over 
the vector of flow counts n = [nj]jgjr G and the vector of queue lengths q = [(7e]eef ^ by 

The following lemma, whose proof is found in Section 5.2, provides a central argument towards 
establishing the strong stability of the fluid model. This lemma will also be of use later in establishing 
the characterization and attractiveness of invariant manifold under critical loading. 

Lemma 3. Let be, respectively, the flow count process and the queue length process of 

a fluid model solution in the set FMS"(T). // L(p) < 1, then for every regular point t S (0,T), 

^L,(n(i),(/(t)) <0. 

Suppose further that L(p) < 1. Then, there exist 5* > and T^, > such that, for all T > T^,, if 
the initial conditions (n(0), ^(0)) satisfy 

L^{niO),q{0)) =1, 

then 

La{n{t),q{t)) < 1 - (5,, for all t e [T^T]. 
5.1. Proof of Theorem 2 

The following lemma provides a sufficient condition for positive recurrence: 

Lemma 4. [21, Theorem 8.13] Let X{-) he an irreducible, aperiodic jump Markov process on a 
countable state space X. Suppose there exists a function F: X — > M+, constants A and e > 0, and 
an integrable stopping time r > such that, for all x with V{x) > A, 

E [V{X{t)) I X{0) = x]< V{x) - eE[r | X{0) = x]. 

If the set {x € X : V{x) < A} is finite and E[V{X(1)) \ A'(O) = x] < oo for all x £ X, then the 
process X{-) is positive recurrent and ergodic. 

Our proof of Theorem 2 relies on establishing the sufficient condition for positive recurrence given 
by Lemma 4, using the stability of the fluid model (Lemma 3). To this end, note that under the 
MWUM-a policy, X{t) = {N{t), Q{t)) eZ^xZ^isa jump Markov process. We shall use 'normed' 
version of the Lyapunov function Lq, defined as 

ein,q)^{L^{n,q))Th , 
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for all (n, q) G M.-^ x M^. The role of V in Lemma 4 will be played by i. Lemma 14 implies that 
there exists constants < Ci < C2 < oo so that, for all (n, q), 

Ci\\{n,q)\\^ < e{n,q) < C2\\(,n,q)\\^. (10) 

Further, for any k > and (n, g), we have that £{Kn, nq) = Ki{n, q). 
Now, consider any sequence of initial states 



such that II 

■^ "lloo — ^ 00 as k — y 00. For each k, consider a system starting at the initial state x . 
Denote the state of the kth system at time t > by 

X\t) ^ {N^{t),Q^it)). 

For each A: G N, define the scaling factor = £(^x^^^'^, and notice that — ?> 00. Fix a time 
horizon T > 0, and for the kth system, consider the scaled state process, defined for t G [0, T] as 

A l^k^rkt) = l{N''{nt),Q''irkt)), 

The descriptor of the scaled system is given, for t G [0,T], by Z^'^''\t) from (5). Let be its 
distribution on D[0,r]. 

From (10), we have that, for any k, 

||;f('^^)(0)||oo = < i/ci. (11) 

Since the set of scaled initial conditions is compact, there must exist a limit point and a convergence 
subsequence. Along this subsequence, by the analysis of Theorem 1 (in particular. Lemma 17) the 
measures {n^^''^} are tight, and therefore, there exists a measure that is a limit point. By 

restricting to a further subsequence, we can assume, without loss of generality, that 
as /c — ?• 00. That is, 

{^^"'Ht))te[o,T] ^ (3W)tG[o,T]' as /c ^ 00, 
with 3(-) G FMS"(T) satisfying the fluid model equations. 

Given a fluid model solution 3(-) of the form (6), denote by (n(-), q{-)) the flow count and queue 
length components. Note that i(^X^^''\0)^ = 1, for all k. By the continuity of £, we have that 
£[n{0), q{0)) = 1; that is La{n{0),q{0)) = 1. Then, from Lemma 3, there exist (5* > and T^, > 
so that, for sufficiently large T, 

£{n{t),q{t)) < 1 - 5*, for aU t G [T„T]. (12) 

Define the functional F : D[0,T] ^ M+ by 

^ {M)tmT]) = f\<tu{t)) dt. 

Since F is a continuous, it follows that 

P ((^^^'=H*))*e[o,r]) ^ P ((3(i))ie[o,T]) ' as A; ^ 00. (13) 
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Further, from (11), and using the fact that the arrival processes are Poisson and the boundedness 
of the rate ahocation pohcies, it follows that for all t € [0, T] and all A; € N, 



{N^''^\t),Q^''^\t)) 



< C.T, 



(14) 



for some constant C3 > 0. From this, the uniform integrability of -F(^^'"'=^(-)) follows. Subsequently, 
from (12) and (13), it follows that 



lim E 



F((z(^-'=)(t))^ 



< l-<5*. 



'te[o,r] 

Equivalently, in terms of the unsealed state process A'(-), we have that 

f-T 



1 



lim , 

fc-S-oo i{x'^) 



T 



;f (0) = 



<l-6,. 



(15) 



To complete the proof of Theorem 2, define C/ to be a random variable that is uniformly distributed 
over [T^:,T]. Define the stopping time r = ^(^(0))C/. Note that 

E[r I X{0) = x] = \l{x){T + T^). 

Then, (15) implies that, for all initial states x € Z:^ x with l{x) sufficiently large, 

E [1{X{t)) I X{<d) = x]< l{x) - eE [r I X{Q) = x] , 

where the constant e is chosen so that < e < 25* /(T + T^,). This satisfies the conditions of 
Lemma 4, and hence completes the proof of Theorem 2. 



5.2. Proof of Lemma 3 

Suppose t S (0,T) is a regular point. We will start by establishing the first part of Lemma 3: if 
L(p) < 1, then ■^La{n{t),q{t)) < 0. To this end, note that 



d_ 



(n°'(t) \ 
^^nf{t)+^q^{tm) • 
^^fP.f ees J 



(16) 



A. 



We consider the terms A„ and Aq separately. 

First, consider the term A.„. For each flow type /, we wish to show that 



nf{t)<qlf^{t){pf-Xf{t)). 



(17) 



By (F2)-(F4), hf{t) = Vf — HfXf{t). There are two cases. If nf{t) = 0, then by (F9), we have that 
Xf{t) = pf, thus both sides of (17) are 0. If n/(t) > 0, then 



(18) 
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Here, the inequality follows from the fact that the function g{z) = z^~'^/{l — a) is concave. For a 
concave function g, g'{y){y — x)< g{y) — g{x), here we have x = Xf{t) and y = pj. Now, since Xf(t) 
is optimal for the rate allocation problem of (F9) and p/ is feasible, we have that 

Combining (18) and (19), we have established (17). Then, we can sum (17) over all / € J-", to obtain 

A„< [r(/,-x(i))]^g"(t). (20) 
Now, consider the term Ag in (16). (F5), (F6), and (F12) imply that 

q{t) = -{I - R^)Us{t)+rx{t). (21) 

From (F7) and (FIO), 

(/ - R^)Us{t)\ ^ g"(t) = V s^(t)7r^(I - i?)g°(t) = max (T^(/ - R)q'^{t). (22) 

Together (21) and (22), imply that 

< [Tx{t)]'^ q^it) - max cj^(I - R)q'^{t). (23) 

Now, combining (16), (20), and (23), we have that 



;La{n{t), q{t)) < (1 + a) ([Fp]^ q^{t) - max a"^(/ - i?)(z"(i)) ■ 

By the definition of \-{p) — PRIMAL(HFp), there exists some s G with l^s = \-{p) and 

Fp < (/ - R^)Us = ^s^{l- i?^)7r. 

This implies that 



(24) 



Hence, by (24), 



[Fp]' g"(t)<L(p)max a' {I - R)q-{t), 



-^L„(n(t),g(t)) <-(l + a)(l-L(p))max a^(/ - i?)g"(t). (25) 

Uit ctSo 



In order to bound the right hand side of (25), we will argue that 



l"rg"(t) 



max cj^(/ - R)q''it) > ^ ' . (26) 



To see this, consider a queue cq ^ £ with maximal length at time t, i.e., Qeoit) = maxg Qeit)- Define 
ei, . . . , ej G £" to be a sequence of distinct queues so that, for each < i < J, packets departing 
from queue go to queue ej+i, and packets departing from queue ej exist the network. By our 
assumption of acyclic routing, such a sequence exists, and since each queue is distinct, J + 1 < \£\. 
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Consider the distinct schedules ttq, . . . , vtj € 5, where, for < i < J, the schedule tTj serves exactly 
the queue Cj — such schedules exist by the monotonicity assumption on the scheduling constraints. 
Clearly, these schedules have weights given by 



(t) if i = J. 



nJ{I-R)q'-{t) 
Averaging over the J + 1 schedules, 

-l—y7rJ{I-R)q''{t) ^ ^. 

i=0 

Since at least one schedule must have a weight that exceeds this average, 

T/T m ^ ^°o(*) ^ maxeg£:g°(t) 
max a [I — R)q [t) > — > — 

(tSiS J + 1 \c\ 

Since 



maxg^(t) > , 

est \C\ 

(26) follows. 

Combining (25) and (26), we obtain, when \-{p) < 1, 

-L4n{t),q{t)) < -(l + a)(l - Up))^^ < 0- (27) 
This establishes the first part of Lemma 3. 

To prove the second part of Lemma 3, we will consider two separate cases over initial conditions 
(n(0),g(0)) with L,(n(0),g(0)) =1: 

(i) l"^gi+"(0) > ei. 

(ii) l"^gi+°(0) < ei. 

Here, ei > is a constant that will be determined shortly. 

For case (i), from the norm inequality in part (i) of Lemma 14, we have that 

< hm\i+a < lk(o)lU, 

Thus, l^q°'{0) > £2 — ^i^" ■ Due to (Fl), q{-) is uniformly Lipschitz continuous, there exists Ti > 
such that l'^g"(t) > £2/2 for ah t € [0,Ti]. From (27), since L(p) < l,we have that 

j^L„{n{t),qit)) < -(1 + a)(l - L(p))^ < 0, 

for all regular t G (0,Ti]. Therefore, there exists 5i > so that 

L„(n(t),(?(t)) < (28) 

for all t>Ti. 
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For case (ii), the argument is more complicated. The basic insight is as follows: if q^^'^{Q) is 
small, then Ag in (16) must be small as well. Further, a good fraction of the flows will be allocated 
the maximal rate (i.e., C), therefore A„ in (16) will be significantly negative. This will leads to 
strictly negative drift in L. We will now formalize this intuition and make an appropriate choice of 
ei > along the way. 

First, by the Lipschitz continuity of g(-), there exists some ri > such that, for all t G [0, ti], 
l^gi+"(t) < 2ei. Note that from (23), the fact that G 5, and x{t) < CI, for ah regular t G (0, n], 

Ag < [rx(t)]^g"(t) < Cl^g"(t). 

By Jensen's inequality, 

l + a 

^l'''Z"(i)j " <^l"'g'+"(0- 
Therefore, for regular t G (0, n], 

Ag < C\T\^ {2ei)^ . (29) 

Again using the fact that l'^q^+"{t) < 2ei for t G [0,ri], it follows that ge(t) < £3 = (2ei)i/i+", 
for all e G f^. Now, consider the set 

J-'^{/G-F : n^(0) >4C£3}. 

By Lipschitz continuity, there exists some < T2 < ri such that, for all / G J-"' and all t G [0,T2], 
nf{t)> 2Ces. Then, by (F9), Xf{t) = C for regular t G (0, T2]. Thus, if / G J", we have, for regular 
t G (0,T2], 

rfPf f^fPf P'fPf 

where we define 

Finally, if / ^ T', we have, for regular t G (0, T2], 

<(t) <(t), , (2C)°p^r"ep^ ^ 

/t^^/W = ^/ - ^ ? ^ /32(2C)"£r , (31) 

PfPf PfPf Pf 

where we define 

P2 = max — . 

Using (16) and (29)-(31), it follows that, for ah t G [0,r2], 



jL^[n{t),q{t)) < (l + a) ( [c\F\^ + h\mcr] 4^" " A E ^/(*) | ' (^2) 



From (27), for all t > 0, La{n{t),q{t)) < 1. For t G [0,r2], l^g(t)^+" < 2ei. Suppose that £1 < 1/4. 
Then, 



™^^/e-^7I7W 4max/g^ ^ 



/eJ-' ■'^■^ m/p5^ -'^-^ A^/py 
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for t € [0, T2]. From part (i) of Lemma 14, 
Then, for regular t G (0,T2], 

j^Lo,{n{t),q{t)) <-(l + a)/3i/33/2<0. 
It follows that there exists ^2 > such that, for all t >T2, 

La{n{t),q{t)) <l- 52. (33) 
Lemma 3 follows from (28) and (33) with J* = min{5i,52} and T^, = max{Ti,T2}. 



6. Critical Loading 

We have established the throughput optimality of the system under the MWUM-a control policy, for 
any a G (0, 00) \ {1}. Thus, this entire family of policies possesses good first order characteristics. 
Further, there may be many other throughput optimal policies outside the class of MWUM-a 
policies. This naturally raises the question of whether there is a 'best' choice of a, and how the 
resulting MWUM-a policy might compare against the universe of all other policies. 

In order to answer these questions, we desire a more refined analysis of policy performance than 
throughput optimality. One way to obtain such an analysis is via the study of a critically loaded 
system, i.e., a system with critically admissible arrival rates. Under a critical loading, fluid model 
solutions take non-trivial values over entire horizon. In contrast, for strictly admissible systems 
under throughput optimal policies, all fluid trajectories go to (cf. Lemma 3). We will employ the 
study of the fluid model solutions of critically loaded systems as a tool for the comparative analysis 
of network control policies. 

In particular, given a vector of flow counts, n € , and the vector of queue lengths, q € M^, 
define the linear cost function 

c(n, q)A^!^ + ^q^ = l^[T diag(M-^)n + q] . (34) 

This cost function is analogous to a 'minimum delay' objective in a packet-level queueing network: 
a cost is incurred for each queued packet, and a cost is also incurred for each outstanding flow, 
proportional to the number of packets that it is expected to generate. 

In this section, we establish fundamental lower bounds that apply to the cost incurred in a 
critically loaded fluid model under any scheduling policy. In Sections 7 and 8, we will compare 
these with the costs incurred by MWUM-a control policies. We shall find that as a — t- 0"*", the cost 
induced by the MWUM-a algorithms improves and becomes close to the algorithm independent 
lower bound 
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6.1. Virtual Resources and Workload 



We begin with some definitions. First, consider the dual of the LP PRIMAL(A) 



DUAL(A) = maximize \^ Q 

subject to n^C < 1 



Note that there is no duality gap, thus the value of PRIMAL(A) is equal to the value of DUAL(A). 

Definition 3 (Virtual Resource). We will call any feasible solution C, G of dual optimization 
problem DUAL(Sr/9) a virtual resource. Suppose the system is critically loaded, i.e., the offered 
load vector p satisfies 



Then, we call a virtual resource that is an optimal solution of DUAL(Hrp) a critical virtual 
resource. 

For a critically loaded system with offered load vector p, let CR(/9) denote the set of all critical 
virtual resources. Note that CR{p) is a bounded polytope and hence possesses finitely many extreme 
points. Let CR*{p) denote the set of extreme points of Q.R{p). 

The following definition captures the amount of 'work' associated with a critical resource, as a 
function of the current state of the system. 

Definition 4 (Workload). Consider a critically loaded system with an offered load vector p and 
a critical virtual resource C, € CR(/9). If the flow count and queue length vectors are given by {n,q), 
the workload associated with the resource ( is defined to be 



6.2. A Lower Bound on Fluid Trajectories 

Consider a critically loaded system with offered load vector p. We claim that the following funda- 
mental lower bound holds on the fluid trajectory under any algorithm. This bound can be thought 
of minimal work-conservation requirement. 

Lemma 5. Consider the fluid model trajectory of system under any scheduling and rate allocation 
policy, with flow count and queue length processes given by (n(-), (7(-)) . Then, for any time t > 
and any critical virtual resource ( G CR{p), 



Proof. Given a time interval [0, T], for any T > 0, consider the fluid model trajectory 3(-) of the 
form (6). By Theorem 1, this fluid trajectory must satisfy the algorithm independent fluid model 
equations (F1)-(F8). By (Fl), the trajectory is Lipschitz continuous and differentiable for almost 
all t G (0,r). For any such regular point t, by (F2)-(F4), we have 



l{p) = PRIMAL(Hr/>) = DUAL(Hrp) = L 



W({n, q) = C~^E [q + r diag(/_f) ^n] . 




(35) 



h{t) 



u — diag{fj.)x{t). 
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Thus, 

rdiag(//-i)n(t) = Tp- Tx{t). (36) 

From (F5)-(F6), we obtain 

q{t) = (l- R^) z{t) - - + T3:{t). (37) 

Adding (36) and (37), we obtain 

q{t) + r diag(/x"^)n(t) = Tp+[l- R^^ (i(t) - Us{t)) . 

Now, multiplying both sides by H = (/ — i?^) ^, we obtain 

E [q{t) + r diag(/i-^)n(i)] = HFp + z{t) - Us{t). (38) 

Consider a critical virtual resource C £ CR(/9). Since C is DUAL(Hrp) optimal, CETp = 1. Taking 
an inner product of (38) with we obtain 

C^E [q{t) + F diag(/i-i)n(t)] = 1 + ^ - C'^Us{t)- (39) 

Now, by (F8), z{-) is non-decreasing. Since C is non-negative, then ('^ z{t) > 0. By (F8), s{t) is 
non-negative. Since C is DUAL(HFp) feasible and from (F7), it follows that C^ns(t) < l'^s(t) = 1. 
Applying these observations to (39), it follows that 

j^w^{n{t),q{t)) = C^H [g(t) +Fdiag(/i-i)n(t)] > 0. 

Given that [n{-),q{-)^ is Lipschitz continuous, the desired result follows immediately. ■ 

Lemma 5 guarantees the conservation of workload under any policy. This motivates the effective 
cost of a state (n, q) € M;^ x M^, defined by the linear program 

c*(n, g) = minimize c{n',q') 

n' ,q' 

subject to Wi;{n',q') >W(;{n,q), V C G CR*(p), (40) 
n eR^, qe . 

The effective cost is the lowest cost of any state with at least as much workload as (n, q). We have 
the following lower bound on the cost achieved under any fluid trajectory: 

Theorem 6. Consider fluid model trajectory of system under any scheduling and rate allocation 
policy, with flow count and queue length processes given by (n(-),g(-)). Then, for any time t > 0, 
the instantaneous cost c(n(t),g(t)) is bounded below according to 

c*{n{0),q{0)) <c{n{t),q{t)). (41) 

Proof. By Lemma 5, if the initial condition of a fluid trajectory satisfies (^^(0), q{0)) = (n, q), then 
(n(t),g(f)) is feasible for (40) for every t>0. The result immediately follows. ■ 
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7. Balanced Systems 

In this section, we will develop a bound on the cost achieved in a fluid model solution under the 
MWUM-a policy. In particular, we will establish that this cost, at any instant of time, is within a 
constant factor of the cost achievable under any policy. The constant factor is uniform across the 
entire fluid trajectory, and relates to a notion of 'balance' of the critical resources of the network 
that we will describe shortly. 

We begin with a preliminary lemma, that provides an upper bound on the cost under the MWUM- 
a policy. This upper bound is closely related to the Lyapunov function introduced earlier for 
studying the system stability. 

Lemma 7. Consider a fluid model trajectory of system under the MWUM-a policy, and denote the 
flow count and queue length processes by (ri(-), . Suppose that the offered load vector p satisfies 
L(/j) < 1. Then, at any time t > 0, it must be that 

c{nit),q{t)) < (l + /3(a))c(n(0),(?(0)), (42) 

where /3(a) — >■ as a — )• 0"'". 

Proof. Recall the Lyapunov function from (9). It follows from Lemma 3 that, so long as L(p) < 1, 

L,(n(t),g(t)) <L«(n(0),g(0)). (43) 
Applying Lemma 14, with p = 1 + a and d = \£\ + \ J-\, 



Now, as a 0+, (ii+'i — ;> 1. Also, 



where z^* = minjz/j and u* = maxjz/j. Thus, as a — > 0^, 1/z^/ 1 uniformly over /. The result 
then follows. ■ 

The following definition is central to our performance guarantee: 

Definition 5 (Balance Factor). Given a system that is critically loaded with offered load vector 
p, define the balance factor as the value of the optimization problem 

'y(p) — minimize c{n',q') 

n,q,n',q' 

subject to W(^{n',q') > W(^{n,q), V C G CR*(p), 
c{n,q) = 1, 
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It is clear that 7(p) > 0, since n',q' > 0. Since there are feasible solutions with {n,q) = {n',q'), 
it is also true that 7(p) < 1. In order to interpret 7(p), assume for the moment that there is 
only a single critical extreme resource C ^ CR*(/>). If we define v = H'''^', then the constraint that 
wc_{n\q') > wc_{n,q) is equivalent to 

[rdiag(/i)~^n' + q'] > v'^ [rdiag(^)-^n + q] . 

In this case, it is clear that the solution to the LP defining 7(p) is given by 

7(p) = (minfe)/(maxi;e). 

e e 

Hence, 'y{p) is the measure of the degree of 'balance' of the influence of the critical resource C across 
buffers in the network. 

In the more general case (i.e., |CR*(/3)| > 1), define the set V = span : ( G CR*(p)}. It is 

not difficult to see that ^{p) > if and only if, for each queue e € V, there exists some v £ V with 
Ve > 0, i.e., if every queue is influenced by some critical resource. We call networks where 7(p) > 
balanced. In the extreme, if 1 G V, then 7(p) = 1. 

The following is the main theorem of this section. It offers a bound on the cost incurred at any 
instant in time under the MWUM-a policy, relative that incurred under any other policy. This 
bound is a function of the balance factor. 

Theorem 8. Consider fluid model trajectory of a critically loaded system under the MWUM-a 
policy and denote the flow count and queue length processes by (n(-),g(-)). Suppose that ^{p) > 0. 
Let {n'{-), q'{-)) be the flow count and queue length policies under an arbitrary policy given the same 
initial conditions, i.e., n(0) = n'(0) and q{0) = q'{0). Then, at any time t >0, it must be that 

c{n{t),qit)) < i±^c(n'(t),g'(t)), (45) 

where /3(a) — > as a 0^ . 

Proof. First, note that if [n(0), q{0)) = 0, i.e., the system is empty, then this holds for alH > (cf. 
Theorem 10). In this case, (45) is immediate. Otherwise, flx t > 0, and set c = c(n(0),g(0)) > 0. 
Define 

in', q') 4 (n'(t), q'{t))/c, (n, q) ^ (n(0), (?(0))/c. 
Using Lemma 5, it is clear that {n,q,n' ,q') is feasible for the LP defining 7(p). Thus, 

c(n(0),g(0)) < -Lc[n'{t),q'{t)). 

The result then follows by applying Lemma 7. ■ 

8. Invariant Manifold 

In Section 7, we proved a constant factor guarantee on the cost of the MWUM-a policy, relative 
to the cost achieved under any other policy. Our bound held point-wise, at every instant of time. 
However, the constant factor of the bound depends on the balance factor, and this could be very 
large. 
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In this section, we consider a different type of analysis. Instead of considering the evolution of the 
fluid model for every time t, we instead examine the asymptotically limiting states of the fluid model 
as t oo. In particular, we characterize these invariant states as fixed points in the solution space 
of an optimization problem. We shall also show that these fixed points are attractive, i.e., starting 
from any initial condition, the fluid trajectory reaches an invariant state. We will quantify time to 
converge to the invariant manifold as a function of the initial conditions of the fluid trajectory. 

This characterization of invariant states is key towards establishing the state space collapse 
property of the system under a heavy traffic limit [4]. Moreover, we shall demonstrate that these 
invariant states are cost optimal as a — > 0"*". In other words, the cost of an invariant state cannot 
be improved by any policy. 

8.1. Optimization Problems 

We start with two useful optimization problems that will be useful in characterizing invariant states 
of the fluid trajectory. We assume that the system is critically loaded. 

Suppose we are given a state (n, q) G M.^ x of, respectively, flow counts and queue lengths. 
Define the optimization problem 

ALGP(n, g) = minimize La{n',q') 

n' ,q' ,t,x,a- 

subject to n' = n + t[i' — diag(/i)x] , 

q' = q + t[Tx- {I - R^) a] , 
n' e R^, q' e M^, t£ M+, 
X € [O,^]-^, aeA. 

Here, recall that A is the scheduling capacity region of the network, defined by (8). Similarly, define 
the optimization problem 

6ALGD(n, g) = minimize La{n\q') 

n' ,q' 

subject to W({n',q') > 'W({n,q), 

VCeCR*(p), 

n' G R^, q' G . 

Intuitively, given a state {n,q), ALGP(n, g) finds a state {n',q') which minimizes the Lyapunov 
function La and can be reached starting from {n,q), using feasible scheduling and rate allocation 
decisions. ALGP(n, g), on the other hand, finds a state {n',q') which minimizes the Lyapunov 
function and has at least as much workload as {n,q). The following result states that ALGP(n,g) 
and 6ALGD(n,q) are equivalent optimization problems: 

Lemma 9. A state {n',q') G M:^ x is feasible for the optimization problem ALGP(n,q) if and 
only if it is feasible for the optimization problem 6ALGD(n, g). 

Proof. First, consider any {n' ,q' ,t,x,a) that is feasible for ALGP(n, g). Note that feasibility for 
ALGP(n,g) implies that 

r diag(^)^^n' + q' >r diag(;u)~^n + q + t 
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Therefore, if C € CR*(p), we have that 

wc{n',q') = wc_{n,q)+t 



Smce (7 S A and Q is feasible for DUAL(Hr/9), we have a < 1. Since ( G CR*(p), we have 
^^HFyO = 1. Therefore, as t > 0, it fohows that 

W(^{n',q') > wc_{n,q). 

That is, {n',q') is 6ALGD(n, g) feasible. 

Next, assume that {n',q') is feasible for bALGD{n,q). Given some t > 0, define 

X ^ diag(^)-i [u - t-\n' - n)] , a = E[Tx- t~\q' - q)] . 

With these definitions, if we establish existence of t > so that < x < CI and o" G A, then 
{n' ,q' ,t, x,cr) is feasible for ALGP{n,q) feasible. 

Note that as i ^ oo, x ^ p. By assumption, < p/ < C, for all f ^ T. Therefore, for t 
sufficiently large, < x < CI. 

Next, we wish to show that, for t sufficiently large, a € A. This requirement is equivalent to 
demonstrating that PRIMAL((t) < 1 and that a > 0. To show that PRIMAL(cr) < 1, note that 
PRIMAL(ct) = DUAL(o-) and suppose that C is feasible for DUAL((t). Then, 

Co = C {EVx-t-\4-q)\ 

= C i^rp - t-'Erdmg{p)-\n' - n) - t'^q' - q)] 
= CETp - t-^ [w^in', q') - w^{n, q)] . 

If C G CR(p), then 

Q^ETp=l, and W(^{n' ,q') — W(^{n,q) 

thus C^o" < 1. On the other hand, if C ^ CR(/?), (^~^ETp < 1. Therefore, in any event, for t sufficiently 
large, DUAL(p) < 1. 

To show that o" > 0, note that 

a = E[Tx-t-\q' -q)] 

= ETp - [Hrdiag(/i)-i(n' - n) + - q)] . 

By assumption, ETp > 0. Therefore, for t sufficiently large enough, o" > 0. ■ 



8.2. Fixed Points: Characterization 

Note that the optimization problem bALGD{n,q) has a convex feasible set with a strictly convex 
and coercive objective function (see, e.g., [1]). By standard arguments from theory of convex opti- 
mization, it follows that an optimal solution exists and is unique. Hence, we can make the following 
definition: 

Definition 6 (Lifting Map). Given a critically scaled system, we define the lifting map A: M'^x 
—7- M.^ X to be the function that maps a state {n, q) to the unique solution of the optimization 
problem 6ALGD(?i, g). 
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The main result of this section is to characterize the invariant states of fluid model as the fixed 
points of lifting map A. 

Theorem 10. A state (n, q) eR^ x is an invariant state of a fluid model solution under the 
MWUM-a policy if and only if it is a fixed point of A, i.e., 

{n,q) = A{n,q). 

Proof. The proof follows by establishing equivalence of the following statements, for every state 
{n,q): 

(i) {n,q) = A{n,q). 

(ii) Any fluid model solution satisfying the initial condition (n(0), q{0)) = (n, q) has (n{t), q{t)^ = 
{n, q) for all t. 

(iii) There exists a fluid model solution with (n{t),q{t)^ = {n,q) for all t. 

(iv) (n, q) satisfy 

(Tp) ^ = max vr^ (/ - i?) , (46) 

TTgS 

PfQiif) =nf, V / G 7". (47) 

(i) (ii): If {n,q) = A{n,q), then it solves bALGD{n,q). Consider a fluid model solution with an 
initial state (n(0),q'(0)) = {n,q). By Lemma 3, it follows that, for all t, La{n{t),q{t)) < La{n,q). 
From the fluid model equations (F1)-(F12), (n(t),g(i)) is ALGP(n,g) feasible, for all t. Therefore, 
it follows that (n(i),g(t)) is an optimal solution of ALGP(n,g), and, by Lemma 9, of 6ALGD(n,g). 
Since 6ALGD(n,g) has {n,q) as its unique solution, it follows that (n(t),g(t)) = {n,q), for all t. 

(ii) (iii): This follows in a straightforward manner by considering the arguments in Theorem 1 
with initial conditions given by {n,q). 

(iii) =^ (iv): Consider a fluid model solution that satisfies (n(t),g(t)) = {n,q), for all t. Then, 
for any regular point t, we have h{t) = and q{t) = 0. Using (F2)-(F4), it follows that x{t) = 
x(t) = p. For any f ^ \i Uf = nf{t) > and Xf{t) = pf < C, then by (F9) it must be that 
Xf{t) = nf{t)/q^(^f^(t). Therefore, PfQi^f) = ^f- Similarly, if nj = 0, it must be that g^^j) = by 
(Fll). 

Now, define H{t) = l'^q^+°'{t). Since q{-) is constant, applying (F5), (F6), (F12), it must be that 
for every regular t, 

= H{t) = q{t)^q''{t) = Wp - (l- R^) Us{t)\ ^ g"(t). 



Applying (F7) and (FIO), 

= (rp)^g" - max 7r^(/ - R)q'^. 



res 



(iv) =^ (i): Suppose {n,q) satisfy (46)-(47). Define {n',q') = A{n,q). Since {n',q') solves the 
optimization problem 6ALGD(?i, g), by Lemma 9, there exists {t,x,a) so that {n' ,q' ,t,x,a) is an 
optimal solution for ALGP(n,g). This solution must satisfy 



n' = n 



+ - diag(/x)x] , q' = q + t 



Fx- I-R^ ]a 
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Now consider the trajectory 

(n(r),g(r)) ^ {n,q) + ^{n' - n,q' - q), Vr G [0,t]. 

Define J to be the Lyapunov function evaluated along this path, i.e., J(t) = La{n{T),q{T)^ . 
Then, 



1 + a 



j(0) nJjuf-flfXf) , g T.. ma 

- > h (rx) ' g'^ - cj ' (/ - 



where S = x — p. 

First, consider Y. Since cr € A, there exists some s G M.'^ with < 1 and a < Us. From the 
monotonicity of S, we can pick s so that a = Us. Therefore, 

cj^(/ - = s^n^(/ - < max 7r^(/ - 

TTgiS 

Then, by (46), it follows that y > 0. Now, consider X, and note that X = by (47) along with 
the fact that 

Thus, we have that J(0) > 0. Since J(t) is a convex function, this implies that J(0) < J{t), 
i.e., La{n, q) < La{n', q'). Due to uniqueness of the optimal solution to 6ALGD(n, q), it follows that 
{n',q') = {n,q). ■ 



8.3. Fixed Points: Attractiveness 



We will now establish the attractiveness of the space of fixed points. Specifically, we will show that 
starting from any initial state, the fiuid trajectory converges (arbitrarily close to) space of fixed 
points, in finite time. 
Given e > 0, define 

Je = {{n,q)eR^xRl : \\{n,q) - A{n,q)\\i < e] . 

In other words, is the set of states (n, q) which are e-approximate fixed points (in an £i-norm 
sense) of the lifting map. Given a fluid trajectory (n(-), q(-)) , define 

hs{n{-),q{-)) ^ inf {t > : {n{s),q{s)) G Je, ^ s > t}. 

In other words, /ie(n(-), g(-)) is the amount of time required for the trajectory (n(-),g'(-)) to reach 
and subsequently remain in the set J',,. 
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Theorem 11. For any e > 0, there exists Hf, > so that if (n(-),g(-)) is a fluid trajectory of the 
MWUM-a policy in a critically loaded system, with initial condition satisfying ||(n(0), g'(0))||oo < 1; 
then 

K{n{-),q{-)) <H,. 
In order to prove Theorem 11, we requhe the following technical lemma: 

Lemma 12. Under the MWUM-a policy, the lifting map A is continuous. Further, A is also 
positively homogeneous, i.e., for all (n, q) € M:^ x and k > 0, 

A.{nn,Kq) = KA{n,q). 

Proof. To establish continuity, it suffices to prove that is {n^,q^) {n,q), then A{n'',q^) — > 
A{7i,q). By definition of A, we have that La{A[n,q)^ < La{n,q). Therefore, since the convergent 
sequence {{n^,q'')} lies in a compact set, the sequence {A{n^,q^)} is contained in a compact set 
also. Define = A{n^,q^) and x = A{n,q). Then, there exists a convergence subsequence of {x^} 
that converges to some x. In what follows, that x = x, and thereby complete the proof of continuity 
of A. 

Suppose that x ^ x. By passing to a subsequence, without loss of generality, assume that 
x'^ X. Since x'^ is feasible for 6ALGD(n'^, g*^), and {n^, q^) — )> (n, q), it follows that x is feasible for 
6ALGD(n, g). Since x is the unique optimal solution to 6ALGD(n, g), we have that La{x) < La{x). 

Now, define 

fc A / W(^{n^,q^) - W(^{n,q) 

e = max 



,CeCR*(p) W(^{1) 

where (x)"*" = max(0, x). Consider x^ = x + e^\. By the definition of e^, it follows that x^ is feasible 
for hlKLQD{n^ ,q^). Then, La{x^) < La{x'^). Now, as fc — ;> oo, La{x^) La{x). Further, — )■ 0, 
so x^ ^ X and La{x^) — )• La{x). Then, La{x) < La{x). By contradiction, this establishes the 
continuity of A. 

The positive homogeneity of A follows directly from the definition of the optimization problem 
6ALGD. ■ 

Proof of Theorem 11. Given 5 > 0, define 

gMJxM^ : L,(n, g) < L,(l)} , 



l^{{n,q)eV 
ls^{{n,q)^V 
lCs^{{n,q)eV 



{n,q) = A(n,g)} , 

\\{n,q)-{n',q')\\i<5, {n',q')el}, 
K{n, q) < K{n\ q'), V (n', q')£V\ Is] 



where K{n, q) = La{n, q) — La[A{n, q)) . The result can be established by showing that the following 
hold: 

(i) K{n{t),q{t)^ is non-increasing in t. 

(ii) For (5 > sufficiently small, X d K-s <ZXs d Je. 

(iii) Starting from any initial condition in T> (this includes all (n, g) with ||(n,(7)||oo < 1), the time 
to hit /C^ is bounded uniformly. 
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In particular, (iii) implies that starting from any state in V, the fluid trajectory hits the set /C5 in 
finite time. By (i), once the trajectory is in set /C^, it remains in that set forever. By (ii), d Je, 
and the result follows. To complete the proof, (i), (ii) and (iii) need to be justified. 

(i) : Since L^^ is a Lyapunov function, then La[n{t),q{t)^ is non-increasing over time under any 
fluid trajectory (cf. Lemma 3). From Lemma 5, the constraints in the optimization problem 
6ALGD(n(t), g(t)) can only become more restrictive over time. Therefore, as the time t increases, the 
cost of the optimal solution of 6ALGD(n(t), g(t)) , i.e., La(A(n(t), q{t))) , can only be non-decreasing. 
Therefore, K{n{t),q{t)^ is non-increasing over time. 

(ii) : First, consider claim X <Z ]Cs for any > 0. A is continuous by Lemma 12. The constraints, 
one for each C, € CR*(/j), in 6ALGD(n,(7) are continuous with respect to {n,q). And (^n{t),q{t)^ 
continuous with over t. Therefore, both functions La{n{t), q(t)) and K (^n{t) , q{t)) are continuous 
with respect to t. Now, D is closed and bounded and Is is open, hence T)\Xs is closed and bounded. 
Therefore, the infimum of the continuous function K is achieved over this set. Since X dX^ for any 
(5 > 0, by the definition of X, this this infimum must be strictly positive. However, over X value of 
K is 0. Therefore, it follows that X C /C^. The claim that Ks C X^ is trivial since if, (n, q) € P \ 
then K(n,q) is greater than or equal to the infimum over that set, hence (n, g) ^ ICs- Finally, to 
establish that Xs C J^e, recall again that A is continuous and, hence, uniformly continuous over D. 
Therefore, for any e > there exists 5 > such that 

\\{n,q) - in',q')\\, < 6 \\A{n,q) - A«g')l|i < e/2. 

Consider any {n,q) G Xs and {n',q') G X with ||(ri,(7) — {n',q')\\i < 6. Then, 

\\in,q) - A{n,q)\\i < \\{n,q) - (n',(/')||i + ||(n',g') - A(n',(?')lli + l|A(n',g') - A(n,g)||i 

< 6 + + e/2 

< e, 

for small enough choice of 5. This completes the proof of (ii). 

(iii) : Here, we shall use Theorem 10. First observe that, by the definition of K and Lemma 5, for 
all < s < i, 

K{n{s),q{s))-K{n{t),q{t)) > L^{n{s),q{s)) - L^{n{t),q{t)) . 

In other words, the decrease in K is at least as much as decrease in L. 

Next, we wish to argue that when the fluid trajectory belongs to the set T> \JCs (i.e., is away 
from the space of fixed points X) then La is strictly decreasing at some minimal rate. To be precise, 
given (n, q) £ 'D \ JCs, suppose that (n(-), g(-)) is a fluid model solution and t a regular point such 
that (n(t),g(t)) = {n,q). We would like to show that 

®(^'^) - TT^ jLa{n{t),q{t)) < -7, (49) 

for some 7 > 0. 
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To this end, for each f £ T, define the function Xf: M.-^ x — > [0, C] as 

Xf{n,q) -- 



Pf if ^/ = 0' 

n'f{t)/q,(f){t) if < n/ < C(?,(j), 

C otherwise. 



Examining (F9), it is clear that this function determines the rate allocated to / at time t, i.e., 
Xf{t) = Xf{n,q). Recalling (16)~(23), we have that 



PfPf 



E 4 (P/ - ^/('^' Q)) + 9)V - ™ - R)q- 

(50) 



/ y2 \ 

E hi - ^/(^' + <f) (^/(^' + [r^r 9" - max - R)q'' 



/GJ- 



Here, for convenience, we define the function T: ^> M by 

T{q) ^ [rp]^ - max {I - R)q'^. 

Now, recall that pf < C for all / € J^. Therefore, it follows that, for all / G J^, if nj > 0, 

Xf{n,q) < Pf Uf/pf < q,(f)it). 

Therefore, for all f £ J-, 

(^"'^"(/)) <0' (51) 

with the inequality being strict if 7^ PfQL{f) nj > 0. 

Since {n,q) G V \ /C^, it can not be a fixed point. Therefore, by part (iv) of the equivalence 
established in Theorem 10, one of the following two conditions holds: 



. T{q) < 0. 

• T{q) = 0, and there exists some f £ T with Uf ^ Pf^i-(f) with Uf > 0. 

Note that the > requirement of the second case follows from (Fll): iif = would imply that 
q'j(j) = 0, hence if iif ^ Pfqi.(f), it must be that Uf 7^ 0. In either of the above two cases, using 
(51) it is easy to see that the right hand side of (50) is strictly negative. However, this does not 
provide a uniform, strictly negative, bound on the drift Ti.liD were established to be a continuous 
function over set V \ ICs, then such a uniform negative bound would follow as the set V \}Cs is 
closed and bounded. However, closer examination (50) reveals that 2) depends on the function Xf, 
thus is not necessarily continuous at the boundary nf = 0, for any / € J-". 

This difficulty is overcome as follows. We will cover D\ICs by a finite collection of closed and 
bounded sets. On each set, we will obtain a bound on the drift D that is continuous on the set 
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as well as strictly negative. Hence, we will conclude that D is uniformly bounded by a negative 
quantity on 'D\ICs, i.e., that (49) holds. The details are given next. 
To begin, define the function 91: IR+ x M+ ^ [0, C] by 



A 



C otherwise. 



For given a vector b = [bj] € {—1, 0, , define Sb as to be the set of (n, q) €^ V \ ICg such that, 
for each /, 

Pf -^{nf,q,(f)) >nf, if b/ = -1, 

-P/l ^ ifbj=0, 

^("/' Hf)) -Pf^ if = 1. 

Clearly Sb is a closed and bounded set, since 9^ is continuous. Given f & J- with bj 7^ 0, define 
the function gf-. Sb ^ by 



9f(ri,q) 



'pJ-^Pf-njT ifb/ = -l, 
{uf + pff-pj ifb/ = l. 

It is easy to check that for all (n, q) € Sb, 

|$H"(n/,g,(;)) -p^l >5/(n,g) >0. (52) 

Further, gf is continuous and gf{n,q) = if and only if Uf = 0. Finally, define the function 
Fb : Sb ^ M as 

FUn,q)^T{q)+ M - g^c/) (/^/ - ^/(n, ?)) 




n\+"gf{n,q)\ 

{c-Pf){pr-c-)n-f, ' y 

We make the following claims: 

(a) Fb is a continuous function over Sb. 

(b) For any (n, g) G Sb, D(n,g) < Fb(n,g). 

(c) For any {n,q) € Sb, Fb(n,g) < 0. 

(a): To establish this claim, it is sufficient to observe that for all / with b^ = 0, x/isa continuous 
function over Sb. To see this, note that if hf = 0, then |9^(n/, — Pf \ < ^^/. Now Xf{n,q) = 
^(^f) <?(,(/)) if > and Xf{n,q) = pj if nf = 0. Therefore, over Sb, we have that Xf{r 
^(^f) *?(,(/)) for all Uf > 0. This establishes continuity of xj for / with b^ = 0. 



n, 



(b): Here, we need to show that for any / with bj 7^ 0, the term in Fb is larger than or equal to 
the corresponding term on the right hand side of (50) in magnitude and preserves the sign. That 
is, 

'< , , / n)+''gf{n,q)\ 

^ - €u) \pf - > -ii^ [iC - Pf){pr - C~nnJ, ^cSf^j ' ^^^^ 
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and 



n. 



n/ 



(54) 



Now, if = then Xf{n,q) = pj and hence the left hand side above of (54) is 0. If nj > 0, 
since hf ^ and thus |9l(nj, gt(/)) — p/| > rij > 0, we have xj ^ pf. Therefore, the left hand side 
of (54) is not equal to 0. Thus, (54) is established. 

To prove (53), we have the following cases: 

• Uf = 0. Here, both sides of (53) are 0. 

• < nj < Cqi^i^fY f^^^^' have 



Pf 



\pf - Xf{n,q)\ 



a 
Mf) 

Pf 



Hf) 



Pf 



\pf-Xf{n,q)\ 



-P/l \pf-H^f^Hf))\ 



^ nfq^^f^gf{n,q) ^ n+'^gf{n,q) 



Pf 

where we have used the fact that 7^ and (52). 
< Cq^f^f) < Uf. Here, since pf < C, we have that 



Pf 



n. 



pf 



Mf) 



\pf - ^fin^1)\ 



q:^fA{c-pf)>{c-pf){py 



(c): Suppose that {n, q) G Sb. We wish to establish that Fb(n, q) < 0. Since (n, q) is not an invariant 
point, by part (iv) of the equivalence established in Theorem 10 and by (Fll), one of the following 
two conditions holds: 

• T{q) < 0. In this case, using (51), clearly fi,{n,q) < 0. 

• T{q) = 0, and there exists some f (z J- with nj 7^ PfQL{f) rij > 0. Here, if bj = 0, then 
since the inequality in (51) must be strict, we have Fb(n, g) < 0. On the other hand, suppose 
that hf ^ 0. Since > 0, we have that 

/ n\+''gf{n,q)\ 
-i- (^(g-P/)(P7°-g->/> '^4. J 

and it follows that Fb(n, q) < 0. 

Now, given claims (a) and (c), it follows that 

sup fh{n,q) < -7b < 0, 
(rt,g)eSb 

for some 7b > 0. Using claim (b), we have that 

S)(n,g) < -7b < 0, 
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for all (n, (?) G Sb- Since P \ /C^ = U 



'be{-i,o,i}- 



lb, we have that, for all {n, q) G V \ ICs 



^{n,q) < -7<0, 



where 



A 

7 = 



be{-i,o,i}-^ 



mm 



7b- 



This completes the proof of Theorem 11. 

■ 

8.4- Fixed Points: Optimality 

The following theorem characterizes the cost associated with an invariant state, relative to the 
effective cost. The effective cost represents the lowest cost achievable under any policy (cf. Theo- 
rem 6). Hence, this result implies that the invariant states of the MWUM-a policy are cost optimal, 
as a — > 0^. 

Theorem 13. Suppose (n*, q*) is an invariant state of a critically loaded system under the MWUM- 
a policy. Then, 



where /3(a) ^0 as a ^ . 

Proof. Suppose (n*, q*) is an invariant state. Define {n', q') to be an optimal solution to the effective 



since (n*, q*) is optimal for 6ALGD(n*, g*), and (n', q') is feasible for 6ALGD(n*, g*). Then, following 
the same argument as in Lemma 7, 



9. Discussion and Future Work 

We have provided a model of a communications network that operates at the packet-level with the 
goal of achieving end-to-end performance at the flow-level. The proposed MWUM-a control policy 
achieves this goal by means of the maximum weight-a packet-level scheduling along with the a-fair 
rate allocation. We established the positive recurrence of the system by means of fluid model when 
the system is underloaded. For the critically loaded fluid model, we established path-wise constant 
factor optimality; the constant factor depends a and the balance factor. 

There are several interesting directions for future work. To start with, by characterizing the 
invariant manifold of the critically loaded fluid model and establishing its attractiveness, the work 
here should lead to the multiplicative state-space collapse property in a relatively straightforward 
manner following the method of Bramson [4]. As the next step, establishing the strong state- 
space collapse property would require bounding the the maximal deviation in the system state 




(55) 



cost LP c*(n*,g*), defined by (40). Clearly 



L^{n\q*)<L^{n',q') 



c{n*,q*) < (1 + /3(a))c(n', q') = (l + /3(a))c*(n*, q*) 



where /3(a) — )■ as a — > 0"*". 



imsart-generic ver. 2009/08/13 file: flow-level-journal.tex date: March 3, 2010 



34 



C. C. Moallemi & D. Shah 



over certain time-horizon. We strongly believe that under MWUM-a control policy for a > 1, this 
should follow from a recently developed Lyapunov function based maximal inequality by Shah, 
Tsitsiklis and Zhong [23]. However, further obtaining a complete characterization of the diffusion 
(heavy traffic) approximation seems to be far more non-trivial question. Finally, the results about 
path-wise constant factor optimality of critically loaded fluid model seem to suggest the possibility 
of such constant factor optimality of MWUM-a control policy under diffusion approximation. 
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Appendix A: Standard Norm Inequalities 

The following lemma provides some standard norm inequalities that are used throughout the paper: 

Lemma 14. Consider a vector a; € M"'. 

(i) IfO<p<q<oo, then \\x\\q < \\x\\p. 
(a) // 1 < p < oo and 1/p + 1/ q = 1, then 

d-^/'i\\x\\i < \\x\\p. 
Appendix B: Justification of the Fluid Model 

In this appendix, we will provide a proof for Theorem 1, which establishes fluid model as the formal 
functional law of large numbers approximation. 

We begin with some technical preliminaries. Fix T > 0. Recall from Section 4.2 that D[0,T] is 
the space of functions from [0,T] to 3, as defined in (6), that are RCLL. This space is equipped 
with the Skorohod metric, defined as 

d(x,y)^inf ||(/.rv||x-(/.oy||, forx,y GD[0,T]. 

Here, ^ is the set of all non-decreasing functions (p: [0, T] [0,T] with (j){0) = and (j){T) = T. 
The norm || • ||° over $ is defined as follows: for ^ € 

4>{t) - Hs) 

t- s 

By i;^ o y refers to the composition y{(j){t)), and for any x € D[0, T], 

||x|| = sup ||x(t)||i, 

te[o,T] 

with II • 111 being the standard ^i-norm over the product space 3. 
For any x G D[0, T] and < s < t < T, define 

'Wx(s,t)= sup ||x(ti) - X(t2)||l. 

S<ti,t2<t 

Further, for any 6 > 0, define 

wU^) = ^ inf max w^{ti^i,ti), 

where S(T, 5) is collection of all 5-sparse decompositions {ti} of [0, T], i.e., Q = t^ < ti < ■ ■ ■ < t^ = 
T with ti — ti-i > 6 for all z > 1. 

It can be easily checked (see [2, Chapter 3]) that the D[0,T] is Polish space under the metric d. 
Let Bt denote the Borel a-algebra on D[0, T] with respect to the topology induced by d. We will be 
interested in probability measures over space (D[0, T],B). We shall utilize the following well-known 
characterization of tightness of measures (see [2, Theorem 13.2]): 
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Theorem 15. The collection of measures {Pg : 9 G G} defined on (D[0, T],^) is tight if and only 
if the following two conditions are satisfied: 

(a) 

lim limsup Pgfx G D[0,r] : ||x|| > ^) = 0. 

(b) For each e > 0, 

lim lim sup Pgfx G D[0,r] : w'^{6) > e) = 0. 
<5^o e&e ^ ' 

We state the following well-known 'concentration' property of Poisson process that shall later 
be useful. It follows from the application of a standard Chernoff bound (see, for example, [9, 
Theorem 2.2.3]). 

Proposition 16. Consider a Poisson process of rate 1. Let N{t) be the number of events of this 
Poisson process in time interval [0,t]. Then, for any 6 G [0,t], 

p(^\N{t)-t\ > <5) < 2exp 

B. 1 . Tightness 

The first step in the proof of Theorem 1 is the following lemma, which establishes tightness of the 
collection of measures associated with the scaled system processes. 

Lemma 17. Under the hypotheses of Theorem 1, for each r > 1, let /i^'') denote the measure of 
G D[0,r]. Then, the collection of measures {/^'-''^ : r > 1} is tight. 

Proof. We will establish tightness by verifying conditions (a) and (b) of Theorem 15. 

First, consider condition (a). It is sufficient to show that for any (5 > 0, there exist K{6) and r{5) 
such that, for K > K{5) and r > r(5), 

/x('')(^x G D[0,r] : M\>K^<5. (56) 
To establish this, fix 5 > 0. By definition, 

l|zM(-)|| = sup (||g«(t)||i + ||zM(t)||i + ||iv(^-)(t)||i + Ws^'Hm, + ||xM(t)||i 

+ ||AM(t)||i + ||I)M(t)||i + pM(t)||i). 

We will bound each component of the system process Z^'^\-). 
First, observe that for any t G [0,T], with probability 1, 

||^(^-)(t)||i + \\S^'\t)\\i + \\X^'\t)\\i < KiT, (57) 

where the constant Ki depends on system dimensions \£\ and \J-\ and on the maximum rate 
allocation C. This is because at most a unit amount of scheduling can be performed per unit time, 
and maximal rate allocated to any flow type is at most C. 

Next, consider the term |[A'^^)(t)||i. For each flow type f G J-, A^J\t) is a Poisson process with 

a time- varying rate that is at most C. Therefore, A^p (T) is bounded above by 1/r times the total 
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number of events of a Poisson process of rate C in time interval [0,rr]. This number of events is 
distributionally equivalent to number of events of Poisson process of rate 1 in interval [0,rrC]. 
Therefore, using Proposition 16, 

P (^5r)(r) > 2Cr) < 2e-^''^^. (58) 
It follows that for r sufficiently large, for any f £ J^, 



P(^4'^(T)>2CTj<^. (59) 
Then, by the union bound, 

P(p(^HT)||i >2|J^|Cr) < A (60) 
Next, consider term ||Q('')(t)||i. Note that 

sup ||QM(t)||i<||QM(0)||i + pM(r)||i. (61) 

te[o,T] 

Now, by hypothesis, we have Q^^\0) — )■ q{0) with probability 1 as r — )■ oo. Therefore, the collection 
of vectors {Q^^\0)} is almost surely bounded, and thus there exists a constant K2 so that, for r 
sufficiently large 

(r) 



i>i^2 <— . (62) 



6_ 
lb 

It follows that for a large enough constant (dependent on K2, C, T, \£\, \ J-\), and for r sufficiently 
large, 

P( sup |IQ(^-)(t)||i>K3) (63) 
\te[o,T] J 

Using very similar arguments (employing Proposition 16), it is possible to bound the Poisson 
processes ||i:>(''')(t)||i and ||AW(t)||i. This, in turn, wih a lead to a bound on \\N^''\t)\\ 1, smce 

E ||ivM(t)||i<||ivM(o)||i + ||AW(r)||i. 

te[o,T] 

Therefore, there exists a constant K4, so that for r sufficiently large, 

P I sup p(^)(t)||i + ||A('-)(t)||i + WN^'Hmi > kA < ^. (64) 



10 



From the discussion above, equations (57), (60), (63), (64), and union bound, it follows that for 
any 6 > 0, there exists constants K(5) and r{6) such that for K > K{6) and r > r{6), we have that 

P ( ll^WfOII >k] <6. 



This completes the verification of (56) or equivalently, condition (a) of Theorem 15. 

Next, consider condition (b) of Theorem 15. For this, it is sufficient to show that for any e > 0, 
there exist 6{e) and r{S{e)) so that for r > r(5(e)), 

P{w'^^r,{6{e))>e) <5{e), 
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with 5{e) — > as e ^ 0. 

To bound w'^^^){5), we need to find an appropriate (5-sparse decomposition {to, ii, • • • , in} of 
[0,T]. For this, we consider a natural decomposition: tj = Z(5, for < i < n, and t„ = T, when 
n = [T/(5J. Then, it follows that 6 < ti — tj_i < 26 for all < i < n. 

We wish to bound — for < s,t < ti, for any < i < n. To this end, first 

note that 

- = ||g('-)(t) - g('')(s)||i + - z^^\s)\\i + ||iv«(t) - iv('-)(s)||i 

+ - 5M(s)||i + ||xM(t) - + ||AW(t) - aM(s)||i (65) 

+ \\D^'\t) - D^''\s)\\i + - A^^\s)\\,. 

As noted earlier, the terms involving S'(^)(-) and X^'^\-) are collectively upper bounded by 

Ki\t — s\, with a system dependent constant i^i, since they are all Lipschitz continuous. Therefore, 
for 5 < £/{10Ki), the sum of these terms in (65) is no more than e/W with probability 1. Next, 
we consider the remaining five terms in (65). As in the justification of condition (a), we will have 
similar argument for all these of five terms. We will focus on the term corresponding to Q^'^^(-)- 
From (3), it follows that 

IIQM(t) _ Q(r){s)\u <\\{i- i?^)n(5W(t) - s«(s)) II, + 1|(/ - i?^)(zW(t) - zW(s)) II, 

+ pM(t)-^M(.)||i. 

Now the first two terms are bounded by K^\t — s\, with the constant dependent on |f |, 
and \ J-\. To see this, note that both S^^\-) and Z^'^\-) are Lipschitz continuous, and the matrices 
(/ — R^)Il and (/ — R^) are finite dimensional (with dimension dependent on \£\, and \J-\), 
and with each entry bounded by constants. It follows that by choosing 5 < e/(10i^5), we have that 

11(1 - i?^)n(sW(t) - 5M(.))||, + ||(/ - i?T)(zW(t) - zW(.))||, < ^, 

with probability 1. Now, to bound the contribution of term ||^('')(t) — ^(^)(s)||i, we can utilize 
arguments used for obtaining (60) with T replaced by |t — s|. Note that here we need to use 
'memory-less' property of the Poisson process crucially. As a conclusion, we obtain that there 
exists a constant Kq so that is (5 < e/{10KQ), then for sufficiently large r, 

P (p«(t) - A«(.)||i < e/io) > 1 - (66) 

From the above discussion, it follows that for any interval as long as we choose 6 

sufficiently small and r sufficiently large, 

P ( sup ligM(t) - gW(.)||i < e/5 ) >l-r^. (67) 

In a very similar manner (using Proposition 16), we obtain the following: there exists a constant 
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K^, so that if (5 = e/Ky and r is sufficiently large, for any < i < n, 



sup 

ti^l<S,t<ti 



|iVM(t)_ArM(,)||^<e/10 >1-^, (68) 



sup ||AM(t) - aM(.)||i < e/io] > 1 - (69) 



sup pM(i) - ^M(.)||i < £/10 ) > 1 - (70) 



sup ||D("^(t) - £'^"^(s)||i < e/10 > 1 - (71) 



lOT 



From (65)-(71) and the discussion above, it follows that there exists a constant K' so that is 
e < K and r is sufficiently large, by a union bound over at most T/5 intervals in the partition, we 
obtain that 

pfmax sup \\Z^''\t)-Z'^''\s)\\i<e\>l-5. 

\ i ti^l<S,t<ti J 

The result follows. ■ 



B.2. Proof of Theorem 1 

We are now ready to prove Theorem 1. 

Given that the collection of measures {//^''•' : r > 1} is tight, for any sequence {r^ : k € N} C M 
with Tfc — )• oo as — > oo, there exists a further subsequence {rk^} and limit point such 
that ^^^'^t' converges weakly to as ^ ^ oo. By restricting to this subsequence, assume that 
these measures are defined on a Polish space, by the Skorohod representation 
theorem, there exists a probability space over which we can define, for all k, random variables 
2('''=)(-) and 3(-) that are distributed according to /u'^''*^^ and respectively, and where the Z^^^\-) 
almost surely converges to 3(-), in the Skorohod metric. We will use this setting to argue that the 
limiting random variable satisfies the appropriate fluid model equations with probability 1. 
Subsequently, we will establish that under an arbitrary control policy, (FMS(T)) = 1, while 
under the MWUM-a policy, (FMS"(r)) = 1. Since jj^'^'i from definition of weak 

convergence, it follows that, under an arbitrary control policy, liminffc_>,oo /^^'^''^'"'^^^(T)) = 1, and 
under the MWUM-a policy, liminffe^^oo /^^'■''HFMS°(r)) = 1, for any e > 0. This wiU imply the 
desired result. 

To this end, we start by establishing that (FMS(r)) = 1. That is, we need to show that 
equations (F1)-(F8) are satisfied. 

(Fl), (F2), (F5): We start with (Fl). Among the components of Z'-''>'\-), Z'-''>^\-), S'^''''\-) and 
X'^'^''\-) are Lipschitz continuous by construction over [0, T]. Since ^(^'=^(-) converges almost surely 
to with respect to the Skorohod metric d, it follows that the corresponding components of 3(-), 
z(-), s(-) and x(-), are Lipschitz continuous. Equivalently, this follows by the Arzela-Ascoli theorem. 

Now, to establish Lipschitz continuity of the other components of 3 we will use the following 
result: 
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Lemma 18. Given a fixed T > 0, consider a Poisson process V with time-varying rate given by 
7(t), for t > 0. Assume there exists constant K > such that j{t) G [0,K] for all t > and that 
7(t) depends only on events that happen up to time t. Consider a sequence {9i} C M+ with 9i oo 
as i oo. Define the scaled process 

v\t) = ^vm, 

for t € [0, T] . Also, define the processes 

y(i) = ^7(M, and f [t) = ffis) ds , 

Jo 

fort G [0,T], where we assume 7'(-) is well-defined. Assume that 'P*(-) converges weakly to V°°{-). 
Then, the sample paths (over [0,T]) ofV°°{-) are Lipschitz continuous with probability 1. Further, 
assume that 7*(-) converges (u.o.c.) to ^{■) over [0,T]. Then, V°°{t) = ^{t), for all t € [0, T]. 

Proof. This is a well-known property of Poisson processes. We describe the key steps of the proof. 
First, using the concentration property of Proposition 16, it can be established that the sample 
paths of scaled Poisson process are approximately Lipschitz. Second, using a variation of the Arzela- 
Ascoli theorem (e.g., see [4, Lemma 4.2] or [33, Lemma 6.3]), it can be established that the limit 
points of such approximately Lipschitz sample paths are in fact Lipschitz. Finally, note the fact 
that number of events under a Poisson process with time- varying rate 7(-) over time interval [0, t] is 
distributionally equivalent to number of events under a unit rate Poisson process over time interval 
[0, Jg 7(s) ds]. As long as 7(-) is uniformly bounded by some constant, say K, the functional strong 
law of large numbers for scaled unit rate Poisson process over [0, KT] can be used to obtain the 
final desired claim. An interested reader may find the details to this argument in many places in 
literature (e.g., [16, Appendix A]). ■ 

Now consider components A^^'='*(-), A^'''^\-), and D^^''''\-). By their construction, these are Poisson 
processes with possibly time-varying rates that are always uniformly bounded. These processes 
converge (over [0, T]) to the corresponding components a(-),a(-) and d{-) of ^{•). Therefore, by 
immediate application of Lemma 18, we obtain that a(-),a(-) and d(-) are Lipschitz continuous. 

Finally, the Lipschitz continuity of n(-) and q{-) is established if (F2) and (F5) hold — this 
is because all of the other components of ^{■) are Lipschitz continuous. To this end, recall that 
equations (2) and (3) are satisfied by scaled system Z^^''^^ for all t G [0, T] by definition. These 
equation are preserved under the almost sure convergence Z^'^''\-) ^{■). Thus, (F2) and (F5) are 
satisfied. 

(F3), (F4), (F6): These equations follow immediately by applying the later part of Lemma 18 for 
Poisson process (possibly time-varying) A(^'fc)(.),A(''fe)(.) and D('^*)(-). 

(F7), (F8): Among remaining equations, first note that (F8) follows because 

Z{'^fc)(.),5(r-fe)(.) and 

X^^f'^-) are non-decreasing and this property is preserved under the almost sure convergence of 
^(^'fe)(.) to A similar argument establishes (F7). 

Now, consider a system that operates under the MWUM-a control policy. We wish to establish 
that /i(°°)(FMS°(T)) = 1. This involves further demonstrating that (F9)-(F12) are satisfied. 
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(F9): Consider any fixed flow type f ^ T and any regular point t G (0,T). Now, if raj(t) = 0, then 
it must be that = 0. This is because of the following argument, utilizing non-negativity of 

nj(-): suppose that either ri/(t) > or ri/(t) < 0. Then, there exist times t~<tort^>t such that 
nf{t~) < or nf{t~^) < — this is a contradiction. Given nf{t) = 0, by (F2) we have a/(t) = df{t). 
By (F3) and (F4), it immediately follows that Xf{t) = Xf(t) = Vf/^f. 

Now, suppose rifit) > 0. Recall that A^|^'°^(-) converges to nf{-) as — )■ oo under the Skorohod 

metric. Therefore, it follows that there exists a 5 > such that, for k sufficiently large, N^j''\s) > 6 
for all s (z [t — 6, t + 6]. We will consider k to be sufficiently large for this to hold. Since t is a regular 
point, Xf{t) is differentiable. Consider any < e < 5. Then, using the fact that Nj^'^\s) > for 
s S [t — 5,t + 6] and the radial invariance property of rate allocation policy, we obtain 

(72) 



x^:''\t + £)-xy''\t) = — Xf{s)ds 

^ ^ rk JrU 




X 



'{Np>) 



argmax — - — (Q^m)"^ I 



x€[0,C] 



a 



Define the function 91: (0, cx)) x M+ ^ [0, C] by 



1 — (y cy 

lH(n, q) = argmax q°^x. 

a:G[0,C] 



1 — a 



It can be easily checked that 

\ n/q if n < Cq, 



C otherwise. 



Therefore, it follows that 9^ is a continuous function. Further, N^^'^\-) and Q[("jj(') continuous 

as functions of time. Therefore, treating d\{^N^j''\-), ^ function of time, it is continuous 

and takes values in [0,C]. Over the bounded interval \t,t + e], it must achieve a minimum and a 
maximum, which we will denote by 9^min(^5 1, s) and 9^max(fc5 respectively. From (72), it follows 
that 

^H„,in(A:, t, e) < ^ < m^^^{k, t, e). (73) 

Now, since (X('''=)(-), A^('"'=)(-), Q('''=)(-)) converges to (x(-), n(-), g(-)) as A; — > oo, it follows (due to 
the appropriate continuity of D'lmiin and IHmax) that 

^min(t, e) < ^f(t + ')^-^f(') < rH^^^(,, (74) 

Here, Dimmit, e) and 5^max(i, s) correspond to the minima and maxima of 5H(n(-), q{-)^ over [t, t + 
Now, taking e — > in (74), invoking the continuity of (n{-),q{-)) and subsequently of 9\, and 
recalling that t is a regular point, we obtain 

x^-"n"{t) 

Xfit) = ^f{i) = argmax — ^ (f^(f){t)x, 

when nf{t) > 0. 
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(FIO): Suppose t € (0, T) is a regular point and consider a schedule vr € iS. Assume there exists a 
schedule o" S 5 with 

7r^(/ - R)q''{t) < a^{I - R)q"{t). (75) 

We wish to establish that ^^^(t) = 0. Since Q^'^'^^-) converges to q{-) as /c — > oo the inequality 
(75) is strict. It follows that there exists 5 > such that for all k sufficiently large and for all 
s G [t-5,t + S\, 

7rT(/ - R)[Q^'-^){s)Y < cj^{I - R){Q^'^\s)Y. 

Then, for unsealed system, the weight of schedule vr is strictly less than the weight of schedule a 
throughout the time-interval [rfc(t — 5),rfc(t + (5)]. Therefore, as per the MWUM-a scheduling policy, 
the schedule vr is never chosen in this time period. That is, for the scaled system, we have 

s^:^\t-6) = s^:'=\t+6). 

Then, as /c — )• cx) 

s^(t + (5) -s^(t-(5) = 0. 

Thus, s^{t) = 0. 

(Fll): Consider a regular point t E (0,T) with nj(t) = 0, for some / G J^. By (F9), we have 
Xf{t) = pf. Suppose that q^(^f){t) > 0. For any ei,e2 > 0, there must exist > so that, if 
s e {t-S,t + d), 

QL{f){s) > ei, and n/(s) < £2- 
Therefore, if k is sufficiently large, we must have 

QJ}J(s) > ei/2, and iv}'''=)(s) < 2e2. 

Equivalently for the unsealed system, 

QL{f){rks) >rkei/2, and A^/(rfcs) < 2rfce2. 

This that, for k sufficiently large and for all s € {t — 5,t -\- 5), the rate allocation in the unsealed 
system must satisfy 

This can be made smaller that Pf/2 by the appropriate choice of £2, and this contradicts the fact 
that Xf{t) = Pf. Therefore, it must be that qi[f){t) = 0. 

(F12): This follows in a straightforward manner from the invariant (for the unsealed system) that 
Z{t) = for aU r G Z+. 
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